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Abstract 

We report effects of local magnetic moment, Gd 3+ , doping (x < 0.3) on superconducting and 
magnetic properties of the closely related Lui-zGd^N^E^C and Yi-zGd^N^E^C series. The su- 
perconducting transition temperature decreases and the heat capacity jump associated with it drops 
rapidly with Gd-doping; qualitative changes with doping are also observed in the temperature- 
dependent upper critical field behavior, and a region of coexistence of superconductivity and spin- 
glass state is delineated on the x — T phase diagram. The evolution of superconducting properties 
can be understood within Abrikosov-Gor'kov theory of magnetic impurities in superconductors 
taking into account the paramagnetic effect on upper critical field with additional contributions 
particular for the family under study. 

PACS numbers: 74.70.Dd, 74.25.Dw, 74.62.Bf, 75.50.Lk 
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I. INTRODUCTION 



With discoveries of new superconducting materials, classical results on effects of non- 
magnetic and magnetic impurities in superconductors^ are being continuously tested and 
augmented. For example, the searches for impurity-induced states in superconductors^ and 
for superconducting quantum critical points^- are few such topics. Unfortunately, in some 
studies of superconductors with magnetic impurities the emphasis is frequently on just the 
superconducting properties whereas the state of magnetic subsystem is often neglected. 

The rare earth-nickel borocarbides (RNi 2 B 2 C, R = rare earth) present a rare opportunity 
to study, within the same family, superconductivity, complex, local moment, magnetism, and 
their coexistence, as well as physics of strongly correlated, heavy fermion materials.- - — In 
this work we concentrate on thermodynamic and magneto-transport properties of LuNi 2 B 2 C 
and YNi 2 B 2 C superconductors with the non-magnetic rare earths (Lu or Y) partially sub- 
stituted by magnetic moment bearing gadolinium. Pure LuNi 2 B 2 C and YNi 2 B 2 C have a 
conveniently high superconducting transition temperatures, T c , and are readily available as 
well characterized single crystals. The details of the superconducting pairing in these ma- 
terials are still debated, with exotic scenarios being examined.— _ — Since the Gd 3+ ion has 
a spherically symmetric, half-filled 4f shell, and therefore virtually no crystal electric field 
effects associated with it, using gadolinium as a magnetic rare earth dopant may simplify the 
problem at hand. Although in resistivity and low field dc magnetic susceptibility the features 
associated with a magnetic subsystem, if located below T c , are often obscured by strong su- 
perconducting signal, it was shown^r— that in this situation heat capacity measurements can 
provide a valuable insight. So far there were several publications, mainly on polycrystalline 
samples, on physical properties of Yi_ :c Gd :r Ni 2 B 2 C^ _ — and Lui_ x Gd z Ni 2 B 2 C^'2i'^i^ solid 
solutions. It is noteworthy that although different studies generally agree on the rate of sup- 
pression of T c (on the pure YNi 2 B 2 C side) and change of the Neel temperature, T N , (on the 
pure GdNi 2 B 2 C side) with x, separation (absence of coexistence) of the superconducting and 



antiferromagnetic order in Yi_ a .Gd a Ni 2 B 2 C near x = 0.3 was alluded to in Ref. 



25 



whereas 



a coexistence of antiferromagnetism and superconductivity at some region of intermediate 



concentrations was suggested in Refs. I20|j23l . Additionally, non-monotonic temperature 
dependence of the upper critical field, if c2 (T), was reported for Luo.88Gd .i2Ni 2 B 2 C.— 
A comparative study of the effects of Gd-doping on T c , H c2 (T) and the state of magnetic 
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sublattice in Lui-^GdrN^E^C and Yi_. E Gd J ;Ni2B2C, has the potential to clarify the effect 
of magnetic impurities on the superconducting state in the rare earth-nickel borocarbides. 

II. EXPERIMENT 

All samples in this study, Lui-zGdzN^E^C and Yi-zGdzN^E^C series, were single crys- 
tals, grown using the Ni 2 B high temperature growth technique^S^. j[ s g rown crystals were 
used for this work. Gd concentrations in both series were evaluated through Curie- Weiss fits 
of the high-temperature part of magnetic susceptibility, that was measured using a Quantum 
Design, Magnetic Property Measurement System (MPMS) SQUID magnetometer. For resis- 
tance measurements a standard, four probe, ac technique (/ = 16 Hz, I = 0.2 — 2 mA) with 
the current flowing in the ab plane, close to I\\a, was used. For these measurements platinum 
wires were attached to the samples using EpoTek H20E silver epoxy and the measurements 
were performed in a Quantum Design, Physical Property Measurement System (PPMS-14) 
instrument with ACT and He-3 options. H c2 (T) data were obtained from temperature- and 
magnetic field-dependent resistance measurements. For these measurements H\\c direction 
of the applied field was kept for all samples. Heat capacity measurements were performed 
in PPMS-14 instrument with He-3 option utilizing the relaxation technique with fitting of 
the whole temperature response of the microcalorimeter. 

III. RESULTS AND DISCUSSION 

A. Heat Capacity and x — T Phase Diagram 

Since the high temperature paramagnetism in the Lui_ rr Gd x Ni 2 B 2 C and Y^^Gd^N^P^C 
series is associated only with a local moment bearing Gd 3+ ion, it was expedient to evaluate 
the real Gd concentration, xcw, by fitting the measured dc susceptibility Xdc = M/H 
(between ~ 150 K and room temperature) with Xdc = xcwC/{T — 0), where 6 is 
the Curie- Weiss temperature, C = (NAplfj)/3k B , Na is the Avogadro number, ks is 
the Boltzmann constant, and p e ff is the effective moment (for Gd 3+ p e ff « 7.94/xg). 
Fig [T] shows experimentally evaluated Gd concentration, xqw as a function of the nom- 
inal concentration, x nom i na i- For Yi_ x Gda;Ni2B2C both concentrations are very close to 
each other (xcw/% nominal = 1-01(1)), whereas the difference is fairly large in the case of 
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Lui_ x .Gd x Ni 2 B 2 C (xcw/ x nominal = 1-31(2)); in both cases the dependence is close to linear 
in the range of concentrations studied. In the rest of the text, the experimentally determined 
Gd concentration will be used. 

Normalized, zero-field, temperature-dependent resistivity data, p(T)/p 300K , for the 
Lui_ x Gd.2;Ni 2 B 2 C and Yi^Gd^N^E^C series are shown in Fig. [2h and [2b. With Gd doping 
the residual resistivity ratio, RRR = p^oml 'Pn, where p n is the normal state resistivity just 
above the superconducting transition, decreases and the superconducting transition tem- 
perature, T c , decreases as well (Fig. [2b, inset). The superconducting critical temperature 
determined from the onset of the resistive superconducting transition for Lui^Gd^N^E^C 
and Yx^Gd^N^E^C is plotted as a function of Gd concentration in Fig. [2b. The T c (x) 
dependence is close to linear with a downturn seen in the case of Y 1 _ z Gd 1 .Ni 2 B 2 C for the 
highest presented doping level. This behavior is consistent with Abrikosov-Gor'kov (AG) 
theory of pairbreaking on magnetic impurities^. The rate of T c suppression is similar for 
two Ri_ x Gd x Ni 2 B 2 C series, being slightly higher for R = Lu. This difference is probably 
due to the additional contribution of the effect of non-magnetic scattering in superconduc- 
tors with anisotropic gap. citemar63a,hoh63a,ope97a,kogl0a Indeed, RRR (that can be, by 
Matthiessen's rule, roughly taken as a caliper of scattering, with lower RRR correspond- 
ing to higher scattering) decreases with Gd doping faster in the case of R = Lu (Fig. [2b, 
inset), that is consistent with larger lattice mismatch (causing stronger scattering) for the 
Gd/Lu (in comparison to Gd/Y) substitution. For comparison, the data for T c (x) evolution 
in Lu(Nii_ x Coa;) 2 B 2 C from Ref. [3J are included in the same plot. It is noteworthy that the 
T c suppression rate is higher for Co-doping to the Ni-site than for Gd-doping to the Lu(Y) 
site, even though among local moment rare earth (e.g. excluding Ce and Yb) Gd 3+ (and 
Eu 2+ ) has the highest de Gennes factor, (gj — 1) 2 J(J + 1), and the strongest T c suppression 
rate.- 1 ^ The reason for such a strong effect of Co-substitution on T c is at least two-fold: 
firstly, Co-substitution for Ni is not isoelectronic, it induces changes in the density of states 
at the Fermi level, therefore causing changes in T c 35M37J]: secondly, for similar concentra- 
tions, x, scattering appears to be stronger for Co-substitution (FigJ2b, inset), thus adding 
to the T c suppressing rate. 

Zero field, temperature-dependent heat capacity, C P (T), was measured for the 
Lui_ x Gd x Ni 2 B 2 C and Yi__ a; Gd :r Ni 2 B 2 C series in order to get additional insight into the evo- 
lution of the magnetic properties with Gd-doping. The results are presented in Fig. El For 



4 



the parent compounds, and several lower Gd concentrations in each series, a jump in C P (T), 
at the superconducting transition temperature is clearly seen. This jump broadens with Gd- 
doping thus the value of AC P at T c was evaluated by the isoentropic construct. Fig. H] shows 
the heat capacity jump inferred from the isoentropic construct for the Lui_ rr Gd x Ni2B 2 C 
and Y 1 _ x Gd x Ni 2 B 2 C series normalized to the value of the jump for the parent compounds, 
LuNi2B 2 C and YM2B2C, respectively, plotted as a function of normalized superconducting 
transition temperature, T c /T c q. As expected, the experimental points lay below the BCS 
law of corresponding states line,— ^ however these points also appear to be below the line 
obtained within the AG theory of pairbreaking from magnetic impurities 4 ^ as well. Simi- 
lar behavior of AC P / AC p0 vs T c /T c0 was observed decades ago for Kondo-impurities (with 
temperature-dependent pair-breaking) in superconductors.— 1 ^ 2 - In our case the dopant, Gd 3+ , 
is a good local magnetic moment ion for which hybridization and Kondo-related physics are 
not expected. There are several possible explanations of such behavior that do not invoke the 
Kondo effect. Qualitatively similar behavior (approximated by AC P oc T 2 ) was observed 
in Yi_ x R x Ni2B2C (R = Gd, Dy, Ho, and Er)^ and was attributed to a combination of 
weak-coupling results of magnetic pairbreaking AG theory with strong coupling corrections. 
Alternatively, a Hartree-Fock approach by Shiba 4 ^ yields a band of possible AC p j AC p q vs 
T c /T c q values that is defined within this approach by the value of the parameter 7, related 
to the strength of spin-flip scattering. For 7—7-1 (weak scattering) the AG results repro- 
duced. The limit of 7 — > describes strong spin-dependent scattering. Our experimental 
data lay close to this 7 — > limit (Fig. H]). Another possible explanation may be a combined 
effect of magnetic and nonmagnetic scattering 4 ^ with a notion that the gap parameter in 
borocarbides is anisotropic. This last possibility is appealing but requires more theoretical 
work due to complexity of the theoretical results and a number of independent parameters 
required for a realistic description. 

Our previous data on the Ybi_ x Gd x Ni 2 B 2 C and Lui_ a; Gd :2: Ni2B2C series^ provide exper- 
imental evidence that for Gd concentration x < 0.3 the long range magnetic order observed 
in pure GdNi 2 B 2 C and the high-Gd end of the series, evolves into a spin glass (SG). A 
broad maximum in heat capacity marked as T max in Fig. [3] is associated with a spin glass 
transition, with T max f» 1.57/, for RKKY spin glasses 4 ^ where Tf is the spin glass freezing 
temperature. 

The other feature in temperature-dependent heat capacity data (Fig. [3]) is a broad min- 
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imum. This minimum exists for all of our x > data and is most probably just a crossover 
between the low temperature magnetism-dominated behavior and high temperature behav- 
ior dominated by electron and phonon contributions. 

Resistivity and heat capacity data together allow us to construct the x — T phase diagram 
for the Lui_ x Gd x Ni 2 B 2 C and Yi^Gd^N^B^C series (Fig. As mentioned above, there 
is a slight difference in T c variation with x between R = Lu and R = Y. The other salient 
temperature, T max , has very similar x - dependence in both cases. It has to be mentioned that 
probing magnetic signatures at temperatures below superconducting transition often is not 
a simple task. In electric/thermoelectric and low field magnetic susceptibility measurements 
the superconducting signal dominates. Magnetic field needed to suppress superconductivity 
might be large enough to alter fragile, low temperature, magnetic state (as it happens e.g. in 
materials with field- induced quantum critical point 47 ), or at a minimum, shift the phase line. 
Zero-field heat capacity measurements clearly reveal (complex) long range magnetic order 
below T c .— In the case of spin-glass transition heat capacity does not have clear anomaly at 
the freezing temperature, Tf, instead a broad maximum is detected at ~ 1.57/. 46 Having this 
in mind, we can approximately outline (by the dotted-dashed line in Fig. |5]) the boundary of 
the spin-glass phase. Since, at least in zero-field resistivity, that was measured in this work 
down to the temperatures below the SG line for several Gd-concentrations, no reentrance 
behavior is observed, superconductivity coexists with the SG state at low temperatures. 
For slightly higher Gd concentrations, after superconductivity is just suppressed, (as it 
was mentioned for Lui_ a: Gd a: Ni 2 B 2 G 1 ^) spin-glass related behavior is observed both in heat 
capacity and magnetic susceptibility. On further Gd-doping, a long range magnetic order is 
established. 

B. Upper Critical Field 

The upper critical field was measured resistively, combining magnetic field-dependent 
data taken at constant temperature and temperature-dependent data taken in fixed magnetic 
field. Examples of such data for Luo.8iGdo.i9Ni 2 B 2 C (i?||c) are shown in Figs. [6] (a) and 
(b). Re-entrant R(T) curves for a few, relatively high, values of magnetic field (Fig. EJ^b)) 
are noteworthy. Results obtained from both data sets are consistent, the resulting H c2 (T) 
curves for two different criteria are shown in Fig. |6]^c). The aforementioned re-entrant R(T) 
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curves are the results of the horizontal (H = constant) cuts through the shallow maximum 
in the H c2 (T). 

The H c2 (T) data for the Lui_ x Gd x Ni 2 B 2 C and Y^^Gd-jN^B^ series are presented in 
Fig. [7J The evolution of the upper critical field behavior with Gd-doping is similar for both 
series: the behavior changes from monotonic with temperature for the parent and lightly- 
doped compounds to the behavior with shallow maximum for higher Gd concentrations. This 
evolution is seen better yet when plotted in normalized coordinates (Fig. [8]). Qualitatively 
similar evolution of H c2 (T) was theoretically described (in dirty limit) by taking into account 
paramagnetic effect.— 1 ^ The use of dirty limit for this materials is consistent with previous 
studies.—. The quantitative description of H c2 (T) within a paramagnetic effect approach 
requires detailed knowledge of the paramagnetic contribution to susceptibility below T c , 
which is a tedious task. On careful examination of Fig. [SJb) we can see that a noticeable 
broad maximum in H c2 (T) is observed for x = 0.14 and x = 0.21, however this maximum 
practically disappears for the next concentration, x = 0.26, for which H c2 (T) is monotonic 
with a tendency to saturation below T/T c « 0.5. For this concentration (at zero field) the T c 
value is close to Tf, the SG freezing temperature (Fig- G2). For spin glasses the paramagnetic 
component of susceptibility decreases below Tf [46j so that H c2 suppression is expected to 
be weaker, in agreement with our observation. Similar arguments were used in Ref. |5l| for 
interpretation of H c2 {T) data below the Neel temperature. 

Fig. [9] presents the slope of the H c2 (T) in the limit of H — > as a function of T c 
in zero field for the Lui-^Gd^N^B^ and Yi^Gd^N^B^C series. The observed behavior 
can be roughly approximated as dH c2 /dT oc T c . It in noteworthy that for several recently 
studied superconductors the behavior is qualitatively different: for Cei_ x La x CoIn 5 dH c2 /dT 
is approximately constant for H\\c and has a factor of two larger absolute value with a slight 
positive slope for i?||a;— for neutron- irradiated MgB 2 dH c2 /dT is approximately independent 
of T c ,— whereas for carbon-doped MgB 2 \dH c2 /dT\ rapidly increases with decrease of T c [54] 
(opposite to what is observed here); for Co-doped LuNi 2 B 2 C the derivative decreases in the 
absolute value only by ~ 20% when T c decreases approximately by half.— 

It is worth mentioning that dH c2 /dT cx T c is predicted for isotropic s-wave materials 
in the clean limit. As discussed above, such description appear not to be pertinent to the 
borocarbides. On the other hand, such proportionality is a property of the AG gapless 
state^^ and is present (at least approximately) in the data from elemental La doped 
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with Gd 



56] (the data in the publication need to be reanalyzed to extract the derivatives). 



Recently, similar behavior was observed in 1111 family of Fe-As superconductors and was 



attributed to pair- 
Following Ref. 



areaking in anisotropic superconductors.— 
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\d(dH c2 / dT) /dT c \ oc n(f)Qk%/h v , where 0o 1S the flux quantum, kg 

0.25 



is the Boltzmann constant and v is the Fermi velocity. Since \d(dH c2 /dT)/dT c \ 
kOe/K 2 (Fig. [9]), the order of magnitude estimate gives v ~ 3 x 10 7 cm/s. This estimate 
is consistent with the values used to describe superconductivity in parent LuNi 2 B 2 C and 
YNi 2 B 2 C.— ^ 



IV. SUMMARY 



Gd-doping of LuNi 2 B 2 C and YNi 2 B 2 C results in T c suppression, consistent with AG 
magnetic pairbreaking with possible additional contribution from non-magnetic scattering 
in materials with anisotropic gaps. For both series T c is suppressed to zero by 30-35% Gd 
substitution. The x — T phase diagram reveals a region of co-existence between super- 
conductivity and a spin-glass state arising from the Gd-magnetism. The evolution of the 
temperature-dependent H C 2 with Gd-doping can be understood by taking into account the 
paramagnetic effect and, for the superconducting sample with highest Gd-concentration in 
this study, Y .74Gd . 2 6Ni 2 B 2 C, by considering temperature dependence of paramagnetic sus- 
ceptibility below the SG freezing temperature. The H c2 derivatives in the limit of H — > are 
approximately linear with zero-field superconducting transition temperatures, in agreement 
with the behavior expected for AG pairbreaking. 

All in all, the Lui_ a: Gd :1 .Ni 2 B 2 C and Yi_ a Gd :1 .Ni 2 B 2 C series present viable systems for 
studies of magnetic pairbreaking in anisotropic superconductors and interplay of supercon- 
ductivity and spin glass state. 
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FIG. 1: (Color online) Gd-concentration evaluated from a Curie- Weiss fit of the high temperature 
susceptibility vs nominal Gd-concentration in Ri-^GdzN^B^C, R = Lu, Y. Dashed lines are linear 
fits with intercept fixed to zero, (see text for details) 
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FIG. 2: (Color online) Normalized resistivity, p/p^ooK for (a) Lui-zGcljjN^E^C (x = 0, 0.055, 
0.13, 0.19, 0.27), and (b) Yi^Gd^rvh^C (x = 0, 0.10, 0.14, 0.21, 0.26, 0.30). Arrows show the 
direction of increasing x, insets: low temperature part of the data. Panel (c): Normalized (to the 
values for the parent compounds) T c as a function of Gd concentration x for Ri_a;Gd x Ni2B2C, R 



= Lu, Y; data for L^Nii-^Coz^E^C from Ref. 
for the same three series. 



34| are included for comparison. Inset: RRR vs x 



13 





FIG. 3: (Color online) Temperature dependent heat capacity for (a) Lui-zGd^N^E^C (x = 0, 
0.055, 0.13, 0.19, 0.27, 0.33), and (b) Yi^Gd^NiaBaC (x = 0, 0.10, 0.14, 0.21, 0.26, 0.30). Arrows 
show examples of how T max and T c are determined. 
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FIG. 4: (Color online) Normalized jump in heat capacity at T c vs normalized T c for the 
Lui-zGd^N^E^C and Yi-^GdzN^E^C series. Dashed lines correspond to BCS law of correspond- 
ing states, Abrikosov-Gor'kov magnetic scattering, and 7 = (strong spin-dependent scattering) 
limit of Shiba's theory.™ See text for more details. 
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FIG. 5: (Color online) T — x phase diagram for the Lui-zGd^N^I^C (filled and partially filled 
symbols) and Yi^Gd^l^E^C (open symbols) series. Symbols: squares - T c from onset and offset 
of the resistive transitions; circles - T c from heat capacity; triangles - T max from heat capacity. 
Dashed lines are guides for the eye. Dotted-dashed line approximates spin-glass phase (crossover) 
line. 
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FIG. 6: (Color online) (a) Examples of magnetic field-dependent resistance of Luo.siGdo.igN^E^C 
single crystal measured at several constant temperatures for H\\c. Onset and offset criteria of 
superconducting transition are illustrated, (b) Examples of temperature-dependent resistance of 
the same sample measured in different applied magnetic fields, (c) Temperature-dependent upper 
critical field of Luo.siGdo.igN^E^C for H\\c. Circles - onset, triangles - offset, open symbols are 
from R(T)\h scans, filled symbols are from R(H)\t scans. 
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FIG. 7: (Color online) Temperature dependent upper critical field (i?||c) for (a) Lui-^Gd^N^EbC 
(x = 0, 0.055, 0.13, 0.19), and (b) Yi-^Gda-NiaBaC (x = 0, 0.10, 0.14, 0.21, 0.26). Circles and 
triangles correspond to the onset and offset criteria respectively. Insets show the same data on a 
semi-log scale. 



18 



CVJ 
O 

X 

CM 
O 

X 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



"> 1 1 r 



^o^9 v o O o 



i 1 1 1 r 

Lu Gd NLB.C 

1-x X 2 2 



H||c 



0.0 



0.0 



> o 
▼ o 



o 
o 



• x=0 

A x=0.055 

x=0.13 
O x=0.19 

i I 



A A 
• A 



o 

T 



0.2 



0.4 0.6 
T/T 



0.8 



1.0 



1.2 



1.0 



0.8 



~i 1 1 r 



~i 1 1 <~ 



Y 1-x Gd x Ni 2 B 2 C 

H lie 



4^ 



* * T 

* 



T 



X 

IN 

1° 



0.6 - 



0.4 - 



0.2 - 



(b) 



0.0 



0.0 



* o 

★ 

▼ 

* 

▼ 

▼ + 



• x=0 

A 0.10 
0.14 
O 0.21 

★ 0.26 

_l i_ 



▼ ★ 

A t 

_ A ★ 

• Air 



0.2 



0.4 0.6 
T/T 



0.8 



1.0 



FIG. 8: (Color online) Normalized to H C 2(T = 
(a) Lui^GdzNiaBaC (x = 0, 0.055, 0.13, 0.19), 
0.26) function of normalized to T C (H = 0) 
criteria are shown. 



0) temperature dependent upper critical field for 
and (b) Yi_ a .Gd a ,Ni 2 B 2 C (x = 0, 0.10, 0.14, 0.21, 
temperature. Data for H\\c obtained using onset 
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FIG. 9: (Color online) dH&jdT in the H — > limit as a function of T c in zero field for the 
Lui-zGdrN^I^C (circles) and Yi-^GdaN^E^C (triangles) series. Dashed line is a guide to the 
eye. 
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